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ABSTRACT 


Initial design study of a low altitude satellite (LASAT) communication sys- 
tem 1s performed. The use of many inexpensive, low-altitude satellites in random 
orbits may be a solution to vulnerabilitv of current military communication sys- 
tems. Statistical study of orbit characteristics 1s performed giving the mean 
number of satellites in view of a ground station and the coverage density of the 
Satellite array. Waveform analysis is performed on a coded and uncoded, 
orthogonal. noncoherent, fast-frequency-hopped M-ary frequency-shift-keved 
signal in a Rician fading channel with optimum partial band jamming. An ana- 
lytical expression for the system probability of bit error is obtained and numerical 
results are generated for various levels of fading. jamming. and diversity. For- 
ward error correction coding 1s applied using convolutional codes and Reed- 


Solomon codes. 
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I. INFRODUCTION 


Current U. S. military communciations are performed using a small number 
of highly sophisticated and expensive satellites in geostationary orbits. The loss 
Of only a few of these satellites will have a serious impact on vital communi- 
cations. Although current technology has not vet prodtced an anti-Satellite 
(ASAT) system that can reach geostationary orbits, such svstems may be devel- 
oped in the future. In addition, ground tracking and control stations are suscep- 
tible to breakdown and attack. To close this vulnerability window many more 
Satellites must be fielded quickly. The problem is the enormous cost of building 
and launching these high-altitude satellites, which is estimated at hundreds of 
millions of dollars each. 

A possible solution to this vulnerabuity in times of crisis is the use of less so- 
phisticated, low-altitude satellites designed to be launched quickly in large num- 
etoe Phese Satellites would be light and inexpensive to build and operate. 
Hundreds of these satellites would be placed in orbit to provide full earth cover- 
age and their sheer numbers would overwhelm any ASAT system. To achieve 
reasonable costs, the satellites will have minimum Station-keeping ability. if any. 
Initial svstem analvsis has indicated that about 250 satellites in random orbits 
flee be needed to provide full earth coverage with inplace reconstitution 
moon 1 he niulti-satellite system (MSS) considered in Ref. | 1s primarly ex- 
pected to be used as a Strategic replacement to the current geostationary system 
and is designed to provide world-wide communications. This configuration 
sera require that the satellites communicate between each other via a packet 
switching svstem. The svstem considered in this thesis is similar to the multi- 
Satellite svstem of Ref. | but would be used for short term tactical communi- 


Savons. | herefore. the satellite crosslinks are unnecessary. 


II. BACKGROUND INFORMATION 


A. LOW-ALTITUDE SATELLITES 

The main purpose of Low-Altitude Satellites (LASAT’s) is to handle crisis 
communications for military commanders. The satellite is expected to operate 
with encrypted data and digitized voice links to local areas, possibly between 
hand-held computers in the field. One advantage of LASAT 1s the small propa- 
gation delay as compared to 250 ms for geostationary Satellites. This will factli- 
tate the speedy transmission of data. and many data protocols designed for 
terrestial radio communications can be modified for LASAT use. 

On the other hand, LASAT coverage 18 very small. and the window of 
communciations will last for only minutes instead of hours. This requires a great 
number of satellites to provide continuous coverage. It 1s envisioned that 
LASAT’s will be deploved randomly in many orbits without station-keeping ca- 
pability and that communications will be through omnidirectional antennas. 
Such a random system requires more Satellites to provide continous coverage than 
does one with well-defined orbits where station keeping is required. Therefore. 
the number of satellites seen by a ground terminal is a random variable, ange. 
ground station signal may be received simultancously by many satellites. Also, 


when used for tactical communications without cross-links, the same ground 
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relayed by many LASAT’s become Satellite multipath signals at the re- 
ceiving ground terminals with near identical power if classical transponders are 
used. In addition. because low-gain antennas will be used by the ground termi- 
nals and LASAT’s, reflection multipath signals are also present. The combina- 


tion of these inultipath effects result m a fading chammer 


B. THE FADING CHANNEL 
Generally, radio waves do not take a direct path from the transmitter to the 
receiver. Terrain. buildings, vegetation. and the ionosphere cause multiple re- 


Hections and refractions of the signal. Therefore. a communicauons channel is. 


tv 


in most cases, a multipath channel with the signals arriving at the recelver from 
Sexaiiety Of directions and with different delays. Furtheremore, the channel is 
Wer, dynamic with the movement of the ionosphere, receiver and transmitters. 
and fohage causing a time variance in the structure of the medium. As a result 
of such time variations, the multipath varies with time. Small changes in path 
length can also result in large phase changes in the signal at the receiver antenna. 

Stein [Ref. 2] describes the system in phasor terms, where the observed re- 
ceived phasor is a vector sum of several phasors with each varing randomly and 
mdividually over a full 2x range. As a result, constructive and destructive inter- 
fememee OCcuis and tlie received signals will have different amplitudes. The fading 
channel, therefore, is characterized by a random time-variant impulse response. 
Thus. the transmission of a pulse will result in a pulse train arriving at the re- 
ceiver with variations in amplitude. If a series of pulses are sent, the pulse train 
as‘ .siated with each transmission will be different as shown in Figure | on page 
ejixel. 3:p. /04]. These time variations are unpredictable to the user and mav 
be considered random which implies that the received signal can be modeled as 
a random process. 

mietire feneral case. the received signal will have both a direct ( or specular) 
component and a number of scattered (or diffuse) components. The scattered 
components arrival times and amplitudes will vary randomly about the mean 
provided by the direct component. As a result. the envelope of the received signal 
Will be similar to the envelope of a non-zero mean signal. This type of channel 
is described as Rician. The amplitude of the received signal has the probability 


@emsiny function (pdr) [Ref. 4: p. 103] 
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figure I. Example of a multipath channel on transmitted pulses: [rom 
Kel 3 pos: 


and all the received energy is via the scattered components, the channel is de- 


scribed as Rayleigh with the pdf of the amplitude given by 





Q Q 

fla) = 5 CA) shoes (Ce Q. (2227) 
G | 26 

The fading characteristic of the channel may be rapid or slow dependent 

upon the relative time and frequency scale of the medium and the signal band- 


width. The coherence bandwidth 1s defined as a iiicasute Ol the siicqueme: 


dependence of the channel. If the channel coherence bandwidth 1s smaller than 
the bandwidth of the transmitted signal, the channel is frequency-selective and 
Mieesienal is Severely distorted. However, if the coherence bandwidth is larger 
than the signaling bandwidth, the channel is frequency-nonselective, and the sig- 
nal is undistorted. A slowly-fading channel 1s one in which the coherence time is 
large and the signal is unchanged over a signaling interval. On the other hand. 
in a rapidly-fading channel, the signal is distorted by changes in amplitude or 
frequency during the signaling interval. 

Unusual conditions in the ionosphere can result in significant scintillations 
fimeartn satellite links. These effects are highly frequency dependent. In mild 
Scintillation the channel 1s typically Rician but at times of strong scintillation 
(deep fading) has been observed to be Rayleigh. The amount of scintillation ts 
@ependent wpon satellite altitude, time of dav, and latitude, with the largest ef- 
fects near the nighttime equatorial zone. Further significant multipath etlects 
that would be involved in a LASAT svstem are dependent upon terrain and the 
Speed of the Satellite and terminals. The relative speed between the Satellite and 
the ground station can be significant with a low earth orbit satellite system. and 
Paiemcrease fading effects. [Ref. 2] 

Mercisiieis one technique that.can be used to overcome the effects of fading. 
With diversity the same information is provided over two or more independently 
fading channels. This lowers the probability that any individual symbol will be 
miected Dy deep fading on each channel. There are a variety of methods of di- 
Meret. cOmbining and thev include the following [Ref. 2]: 


@espaced antennas 
Antenna spacing of ten Wavelengths or greater can result in independent 
channels. This technique. however. is not well suited to satellite UHF com- 
munications because of the long wavelengths. 


Wetec ueiicy diversity 
The use of wideband frequency hopping (fast hopping) on fading chan- 
nels. with the same information carried on different hops, tends to result in 
independent states of fading. 


TAN 


e Multipath diversity 
A signal having a bandwidth greater than the coherence bandwidth of 
the channel will resolve multipath components and thus, provide the the re- 
ceiver with several independent fading signal paths. This is an inherent 
characteristic of direct Sequence spread/spectuinn ac tems. 


> Time divenci 

The same data 1s repetitively transmitted at intervals that exceed tie 
coherence time of the channel which results in independent fading condi- 

lions. 5 
Coding and interleaving are other techniques which are similar to time di- 
versity in providing redundancy. This is particularly true with binary Timea 
block codes. However, convolutional codes and other non-binarv linear codes 
provide satisfactory results. A more detailed coverage of fading channels is pre- 


Scented mokGie-. wands: 


C. FREQUENCY HOPPING SPREAD SPECTRUM! 

The previous section discussed the effects of the fading channel and indicated 
some methods available to overcome the resulting degradation in svstem per- 
formance. Fast frequenev hopping spread spectrum, with the same information 
carried on different hops. 1s one form of frequency diversity that may Dewiseq 
If sequential hops exceed the coherence bandwidth of the channel. then fadinauren 
each hop frequency is independent. The selection of the frequency slots in each 
signaling interval 1s made pseudo-random according to the output of a pseudo- 
noise generator. 

Due to rapid signal phase changes which occur within a fading channel and 
the difficulty in building and using coherent demodulators for frequency shift 
keved systems. most fast frequency hopping systems use noncoherent. orthogonal 
modulation schemes. A common data modulation for FH systems is M-arv fre- 
quency shift keving (MFSK). were A = log, Wf bits are used to determine wines 
of .s/ frequency tones will be uSec for carrier modulation. [he tones are usm 
spaced far enough apart so that the transmitted sicnals are orthosenale ia 
implies that the frequency spacing be at least 1/A 7, Hizesvaere® 7, is thie pie 


ration ane J.= A 7,48 the Curae -ometen «moor 
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When fast frequency hopping is used with a M-ary FSK system, each of the 
Mf symbols is subdivided into Z chips. At each chip, the MFSK modulator out- 
put is hopped to a different frequency. Since the chip duration 7, is shorter than 
the data modulator output symbol! duration 7,. the minimum tone spacing for 
Memmoeonal signals is B=1/7,=L/KT,. The total system bandwidth then be- 
Pommes i = wb = mLjAlj/AT7,, Where m is an arbitary integer that defines the 
hopping bandwidth. : 

A typical frequency hopping system is Shown in Figure 20n page 8. At the 
heeelver, an identical PN generator synchronized with the received signal removes 
the pseudo-random frequency translation by mixing the synthesizer output with 
the received signal [Ref. 5]. For an orthogonal. noncoherent system, demodu- 
lation is usually accomplished through square-law, or envelope detectors. The 
Sevector OULPUL can then make a hard decision at each chip interval as to which 
of the M-channels contains the signal and make an estimate on all £ chip deci- 
Solis. OF 1t Can Sum the energy of 1 chips on each channel and make a soft deci- 
eye based on the iargest value. 

Although the fast frequency hopping svstem 1s effective in a fading channel. 
partial band jamming can seriously degrade the signal when it hops in and out 
Sreeie jamming band. fhe partial band jammer concentrates its power over a 
small fraction 77 of the svstem bandwidth. On some hops, the signal will be at a 
feemency which is jammed, while on other hops, the signal will be at a frequency 
mete ie ONY iiterference is thermal noise. A svstem which can measure the 
amount of thermal and jamming noise power per hop can provide this side in- 
foraraiion to normalize the measurements obtained at the detector outputs 
(chapter IV). Worst case jamming occurs when the jammer optimizes 77 to max1- 
mize the system bit error probability. In svstem design it is normally assumed 
that the ma ines Knowledge Of tjesexsienmestate andus able to optimize 7 . 


Therefore. performance calculations provide an upper bound on error rates. 


D. SUMMARY OF RESEARCH 
This thesis covers the communications analvsis of LASAT svstem design. 


memecd i 1S study, arc orbit characteristics of the system. and the waveform 
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Figure 2.) Typical Frequency hopped spread spectrum system 


analysis of a proposed fast frequency hopped M-ary frequency shift keved mod- 
Ulan slic tne: 
I. Orbit statistics 

Chapter II provides a study of the LASAT orbit distribution and cover- 
age. The first section of the chapter looks into the number of satellites in view 
of a ground station. An analytical approach was used to derive the mean value 
and the standard deviation of the number in view. A computer simulation veri- 
fied the results. The second section of chapter II) discusses the covcrace OMe 


LASAT relative to a ground terminal. The mean coverage area and the standane 


deviation in which connectivity can be established with other ground terminals 
was found. A computer simulation verified the results. 
2. Waveforin selection 

Waveform selection involves spread spectrum techniques, carrier modu- 
lation, error correction coding, and diversity. Chapter IV presents the derivation 
of the bit error probability for uncoded and coded fast-frequency-hopping M-ary 
frequency-shift-keving schemes in a Rician fading channel 4vith optimum partial 
band jamming. Numerical results for various levels of diversity and fading are 
provided for both the coded and uncoded systems. Convolutional codes and 
Reed-Solomon codes were used to provide forward error correction coding. 


Chapter V presents conclusions and recommendations for further research. 
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Ill. LASAT SYSTEM 


A. DISTRIBUTION OF LASAT 

Since the LASAT communication system is envisioned as consisting of a large 
number of low earth orbiting satellites distributed randomly throughout many 
orbits, it is expected that the number of satellites in range of a ground station Will 
also be random. The distribution of satellites within view of a ground station Will 
be developed here both theoretically and numerically. The ground station is as- 
sumed to have an omnidirectional antenna that will reach satellites with an ele- 
vation angle of zero degrees. A general structure of L Satellites randomly 
distri>uted within an orbit altitude range /, to A, will be assumed. 

i. Theoretical development 


The total volume 7 enclosed by the orbital radii A, and A, 1s 


n(R3 — Rj) (3. 


ty | 
{y2 
—" 
ee” 


Ky =RK, +h, 


(J 
ke = 1 + ts 


with A, = 6378.14 km. (3443.93 NM), the radius of the earth. The volume den- 
Sity Of 7 Satellites Viti ti. Olumie ds 

; ; 3 Ale 

p= 15) SSareliivess aie (3.3) 


The velume of the orbit portion seen by a sround station can be derived with tie 
help of Figure 3 on page 12. The volume of the spheniealescctor 12 dei megmer 


points ACBO with radius A, is 
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Eaere @ 15 the radius of the base and a*= RA: — K?. The volume of the segment 
and base ABC 13s the difference between the volumes of the spherical sector and 
Meemieht circular cone. 
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By a similar method, the portion of volume seen by the ground Station below /, 
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Figure 3. Geometry of satellite orbit distribution. 


"y= ZL AjFBR, — hy), (3.8) 


and the subvolume rv of the orbital shell of thickness 4, —/h, seen by the ground 


station Is 
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Therefore, on average, the number of satellites seen by the ground station will be 


Bieeensity Of Satcilites times the subvolume y, or by using (3.3) and (3.9) 
Te — Vil. (3.10) 


Biilesamove result can also be Shown Statistically with the help of the dc- 
vVclopment by Larson [Ref. 6: pp. 143-149]. The distribution of satellites 
throughout the volume J’ is random and uniform and one can assume that a 
Polume v can be divided into 77 non-overlapping subvolumes of size Ar. where 
Ay = v/72 1s Verv small. It is therefore reasonable to assume that there will or will 
not be a satellite in each of these small volumes and that the probability that ex- 
Semeeee@uc Satellite will be in cach subvolume is pAv = pv/7 for cach. where p is 
fi@eecensity Of L Satellites distributed throughout the total volume J. Further- 
iene. il Can be assumed that the presence of a satellite in cach subvolume 1s in- 
dependent of the other subvolumes. Thus the volume ry has been subdivided into 
Peepcated independent Bernoulli trials. each of a volume Av =v//7 with the 
Bie aoity Oi a Satellite being present equal to pAv = pv/m. If mis the number 
of satellites in the subvolume v. then 7 1S approximately a binomial random vari- 


able with probabiltity density function f\(7) where 


Hn , man 


cae a 
=( NG a - | (1 | ey Oe a aed (3.11) 


~ 


= 0. otherwise. 


If the expression is taken to the limit 77 > co (thus the individual subvolumes are 
Shrinking to zero), the exact probability function for the Poisson random variable 


results. Phatis 


: (3.12) 


Therefore, the distribution of the number of satellites within a volume v can be 
represented by the volumetric Poisson distribution with a parameter pr. The 


well-known mean # and variance o- of the Poisson process 1s 


R= pvaat, B48) 
and 
oh= pvt (3.14) 


2. Numerical simulation 
To generate a sample set representative of the distribution over the orbital 
Volume. a uniform random sequence of Z points between O and F is generated. 
To find the number of satellites seen by the ground station for a particular sample 
set of L points between 0 and J’, the number of points between 0 anda 
counted. where I’ and vy the orbital volume and subvolume defined by (3.1) and 
(3.9). When a large number of sample sets are run and the results averageda ae 


mean number of satellites and the variance co: of the samples can be found from 


i=) 2, (3.15) 


and 


Where .\f is the number of sample sets and 7, is the number of satellites found 
Within the subvolume v in the i th sample set. : 

Sample means were found by generating first 2000, and then 5000 sample 
sets of 250 and 150 points untformily distributed as described above. The num- 
ber of satellites found in the sub-area were then averaged and compared against 
Biemmican found by equation (3.13). Table ! shows the numerical results for 
A, = 550 km (300 NM) and #2. = 750 km (400 NM) with 2000 and 5000 samples. 
The numerical results compare favorably with the theoretical results obtained bv 


feoand (3.16). Differences are due to the pseudo-randomness of the uniform 


Mmemioer Sencrator used. 


Table 1. MEAN AND VARIANCE OF SATELLITE DISTRIBUTION. 


[| Theory | 2000 Samples | 5000 Samples 
Me | as | a>. | a 
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oie 


Mehen: aveexperiiment described by a Poisson process is repeated a sreat 
mumteecr '/ Of times. then it can be expected that the number .\/, of times that 


eet 7? Satellites are observed will be (Ref. 7: p. 158] 
VE MSP). ical) 


ieee criiy this. the number of satellites found within the subvolume for the run 
ewer S000 Sample scts was plotted versus the number of times (occurences) A/ 


n 


Meieomeariieuler number of satellites were within view (see Figure 4+ on page 


— 
tsa 


17). This curve is compared to a Poisson distribution with the parameter pv as 
siven by G.13), and shown in Pablo sixsee micire sr omecccesin a: 

Table 2 shows the maximum number of satellites 72, secn by the proume 
station with a .999 and .990 probability. Since the distribution is Poisson, the 
probability P that there will be no greater than n,,, Satellites seen by the ground 


Station 1S 


i=) Sh 
max 


(py)" —pv (3.18) 
7 n! ete 


The probability P, that a ground station sees no Satellites overhead is (3.18) with 


Heed 


jo ot 
JP ce (Sa 


and 1s also given in Table 2 


Table 2. NUNIBER OF SATELLITES SEEN BY GROUND STATION 





The above analvsis shows the expected number of Satellites a ground Sta- 
tion may communicate with. For a tramsmitting:rece;rving pair of ground 
Stations. the number of satellites that may be used for communicating 1s depend- 
cnt upon a joint pdf between the number of satellites, the orbital altitude and the 
distance between the stations. An analvsis of the satellite coverage will be cov- 
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Figure 4. 


ideal Poisson distribution with parameter pv = 11.54. 
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Distribution of the number of satellites within view: 


NUMBER OF SATELLITES = 250 
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Compared with 


By (COVERAGE 
1. Theoretical development 
Now that the mean number of satellites seen by a ground station has been 
established, it is desirable to develope a distribution for the coverage areca of the 
Satellites. A satellite on the horizon of the base Station has a sub-point that is an 


angle « away from the base station (sec Figure 5 on page 19), where 
—| 
t= COSM xe nor (3.20) 


The Satellite at point S can communicate an additional distance of arc length 
R.w to a user terminal at point UT. Therefore. the base station at point BS has 
a maximum communication range of arc length 2¢k, out to the user terminal 
through the satellite. This gives a maximum coverage area defined by the solid 
angle of 4% or 4cos-'(A,/R,) steradians, which for an average altitude of 
h = (hy+h,)/2 = 650 km (350 NM) is an area of 
> ol = 

Re COS RK, aaa 
17.63 x 10° km’. 


sus 
| 


(3.21) 


Figure 6 on page 20 shows an cXample of the coverage area providedmiag 
a number of satellites. The spherical caps of arca covered by the satellites are 
Shown as circles of an equivelent area. In this diagram there are three satellites. 
Si. S2. and S3. uniformly spaced over the total view area of the base siamam 
(Spherical cap of a radius of arc length «&, from the base station). Each satellite 
has a range of arc length aA, also. User terminals, 11, 12, 73. '4) and {Seam 
also randomly spaced over an area out to a range of arc length 2uA, from the base 
Station. One can see from this Figure 6 that the user terminals have differing 
numbers of satellites through which they can communicate to the base station. 
In this example. the user terminals can communicate as follows: T+ through three 
satellites. T2 through iwo satellites. T] and T5 through one satellite cach. ana 
T3 has no Satelite in range. In general, as the ranee Gem cen the user tes ames 


and the base station increases. there are fewer sateliites throwen Ws) icles 


BS = Base station 
UT = User Terminal 
S = Satellite 





Figure 5. Geometry of satellite coverage. 


communicate. If more satellites were added to the system, the communication 
ablilitv would increase. 

iiieweovcraccscensity 10M the totalenumber of satellites L im the orbital 
volume V’ needs to be determined. This density represents the number of satel- 
lites that can communicate with both a user terminal and the base station as a 
function of the distance r between them for a given number n of satellites that are 
Within view of the base station. The satellites are uniformly distributed through 
tie orbital volume and their different altitudes will result in coverage areas of 


Seercmt sizes, This adds another random variable to the process. Since the 
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Figure 6. Example of satellite communication range. 


altitude difference 1s assumed to be small compared with the coverage area of 
each satellite, to simplify calculations without greatly affecting the results, it is 
safe to assume that all the satellites are at the same altitudevand that thesemems 
are circular. The mean altitude /# will be used in calculations. The area covered 
by the satellites can be modeled as an %-Y Cartesian srid with the base Siamen 
at the origin. The positions of the individual satellites are represented bv coor- 
dinate points x, and 3,, which are uniformily distributed from ~ «R, to oR, along 


the X and Y axes. Therefore, the pdfs for the random variables x and y are 


frlx) = 1)20R,, — oR, <x <a, 


Y 
fy) = 1/20, aR Sy SAR, (3.22) 


with a zero mean and a second moment E[x) = Ely) = quae = (a) 3. 
For a given number of satellites 1, there are n sets of randomly distributed 


yil= |2Qy.., cE achmaiaie 


HS, uF 


coordinates describing the satellite positions, x 


satellites covers an area defined by the radius of are length «&, around the 


Satellite subpoint. The average distance from the base station to a satellite 


Subpoint can be described by two new random variables 


n 
x= fi 
i=] 
' (3.23) 
z — ule Y > 
ae Me 
i= 


iemetess, and z, are the average A and Y axes coordinates for the satellites. Since 
nis a Poisson random variable, the distribution of z, and z, will be conditioned 


on#. The conditional probability density functions of z, and z, are given by 


Sf Z,.| yx Jn) = ~ [Py (fx (xa) Ly (%n)] 


_— (3.24) 
fe eyl) =U OWT” he, Wdh 


where the * is the convolution function. The central limit theory provides that if 
#1S Suificiently large (greater than five). the sum of the random variables ap- 
proacn a Gaussian distribution. Thus. z, and z, can be approximated as inde- 
pendent indentically distributed Gaussian random variables with zero mean and 


a variance of o; 














eee ae 
- nad = Ine- ory 2G7 5 
ip Spee heey Z 
3 (3.25) 
fz | qlZy| 0) = = ; ae = 
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where 








2 2 2 
o, = E[z,] = E[z,] 
as ai | nN n 
=E| +) ) xx l=4) > Elsa 
MN j=lj=l 25 = (3.26) 
Nn 
| pe (aR) 
a? > Bi Sn 
i=] 
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Since z, and z, are independent. their joint density function is given by 





n= te inlay | n)tz, Fics: | 72) 
( ze tr ee (3.27) 
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tv 


The base station can communicate with an area in the Shape of a cmee 
with radius 2%A, and the distribution can be represented as a function of the 
distance r of the user terminal to the base station. To do this, define a new ran- 


dom variable 


ae (3.28) 


wihieli 18 the range from the base stattomo tne Use terminal, aac 


Fri ft li) = | J2Z,in'ee% | nde, dz, 
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not 
The region C is a circle of radius r. The above integral can be easily solved using 


polan coordimates [Reivo.p. 7), eet 


z.=rcos @ 
- =rsin 0 (3.30) 


dz, dz, = rdrd€, 


a 


i fom 


then 
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The unconditioned density function fp(r) can be found by multipling by (3.12) 


and integrating with respect to 7 


fr(t) = | . Sri lt laty(ajdn 


7 | 31 a —3rn (pv) pa) 
eee UG) a 


2. Numerical simujation 
A numerical simulation of the coverage problem was run to verifv the 
Meove results. in the model, the average number of satellites that were in view 
of both a user terminal and the base station were found for various ranges from 
the base station. This was done bv using a sampling raster of points spaced 45 
degrees apart On concentric rings around the base station. For each test point, 
tae number of Satellites that were within its view range was counted (a Satellite 


Peesewithin range jf the distance to it from the base station was less than @A,). 


A total was found for all the test points on each ring and the number was divided 
bv the number of test points on the ring (nine) to get the average number of sat- 
ellites that could be seen from that range. 

The sub-satellite points were generated by using a uniform random num- 
ber generator to provide X and Y axis points that were within a radial disiames 
of one from the base station (the value of «AR, was normaljzed to one). The set 
of test points Were generated in concentric rings about the origin out to a radius 
of two. For each test point, the number of sub-satellite points within a radial 
distance of one were counted (since all sub-satellite points were within a radial 
distance of one of the base station it was unnecessary to test the base station to 
Satellite distance). The total number of satellites in view for each test ring was 
then averaged to provide an average number of satellites that can communicate 
with both the base Station and TieMiser terminaiton viiatencmace 

The test results for three different runs of 13 satellites each are plowmeds 
along with the theoretical result from equation (3.30), in Figure 7 on pa 
The theoretical results were multiplied bv the the number of satellites in the array 
to normalize them with the simulation count. The exponential character of the 
numerical density 1s clearly seen and matches very closely to ine theorciedians 
sults. The difference out near the edge of the range 1s due to the approxima 
of the Gaussian function by the central limit theory near the tails of the distrib- 


ution. 
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Figure 7, Coverage density of satellites: As a function of the distance between 


ternunals. 


ho 
ca” 


IV. WAVEFORM ANALYSIS 


A. UNCODED FH/MFSK IN A RICIAN FADING CHANNEL WITH 
PARTIAL BAND JAMMING 

The LASAT communications system under consideration must be able to 
overcome a variety of performance degrading phenomena. Chief among these 
are: Additive White Gaussian Noise (AWGN), jamming and interference, and 
fading. Partial band jamming is the optimal jamming Strategy most likely to be 
encountered. Interference may be a result of, amons other things, othememaa. 
munications svstems. multiple access Strategies. or multiple transmitters in the 
link. For the LASAT communications system the most likely fading encountered 
will be Rician with a direct line-of-sight component and a random. faded com- 
ponent. 

The communications system waveform under consideration in this thesis will 
be required to overcome the above conditions. The system chosen for study ts a 
fast-frequency-hopped M_-ary frequenev-shift-keved (FH, MFSK) system. It 1s 
assumed that the channel will be jammed by an intentional jammer whose power 
resource 1s AWGN. Channel fading is assumed to be Rician fading with a ran- 
don faded component due 10 tenium Teller tiem 

Lindsey [Ref. 9] studied a NIFSK system with diversity in a fading chanme: 
and showed that diversity is effective in overcomming the effects of iading Gnutine 
channel. Lec et al. [Ref. 5. 10] considered a special Case of the Lindsey svstem 
with the diversity apphed by fast hopping the MFSK signal. Jamming and cod- 
Ing Were applied tn the analvsis but fading was not. The study done here 1s a 
continuation of the work by Lindsey and Lee et al., by considering jammung in a 
Rician fading channel of a FH MFSK svstem. 

The transmitter receiver structure similar to Fig. 2 of Ref. 5 is used. A block 
diagram of the transmitter 1s shown in Figure 8 on page 28. The binary dara 


provided at a rate A, = 1/7,.where 7; is the data bit Jenctnsia a led secre 
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Byeer wv iich promdes Af-= 24 symbols at a rate R,= R,/K Hz. There are L hops 
@emesvnibol, Each of the svmbols is then hopped at a rate XR, = LR, Hz. Fre- 
quency hopping 1s emploved over an allowed system bandwidth IV Hz. The 
carrier 18 hopped to a new frequency at least once per svmbol duration. Thus, the 
Bandwidth of the transmitted radio-frequency pulse is 8=2/2, where 
t=17,/L,7,=symbol duration and the bandwidth of the M-ary cluster will be 
fi B= 2 \f/t. 

iimeeceiver iS Shown in Figure 9 on page 29. The dehopped signal r(f) is fed 
into M channels. The output of the square-law detectors are sampled once every 
me seconds 10 produce the quantities {x}, i= 1,2,....../;4 = 1,2,..,L. The 
Serector OUtput Samples are multiplied by the weights 1/e; , which are outputs 
of a noise onlv channel estimator, to produce {z,} and are summed over the index 
felon each i to form the decision statistics {2}. [Ref. 5] 

leit is considered that any portion of the System bandwidth MW’ has an 
equallv-likely possibility of being jammed over the jamming bandwidts FH’/1, 
where 7 1s the jamming fraction, then a single cell of bandwidth B has a proba- 
bility # of being jammed and a probability ]- 7 of not being jammed. Therefore, 


after dehopping. the received signal in the Ath hop can be represented as 


(4.1) 


| fl + {t) + j,(0) With probability 7 
(2) | | 


7 Mele) ee, with probability 1 — 7 


where 7,(t) and j,({t) are the independent thermal noise and jamming noise com- 
moments, respectively. 
Without loss of generality it 1s assumed that the Ath hop signal is present at 


channel ! and is represented by 
Sf) = a, COS(@) tf + 0,). (4.2) 


aid has a transmitted cnerev in the Ath hop given by 
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Figure 9. 





‘Jie ak =e (13) 








Where it is assumed that all transmitted signals have equal energy E. The pa- 
rameter a, in (4.2) is the Rician variable representing the signal amplitude. The 
signal s,(t) has both diffuse and direct (Specular) components. Therefore. as dis- 
cussed in section 2.B the coefficients are Rician distributed and have the density 


function f, (a) given by (2.1) 





=) >) 
Ar Gio or Cate QT Je ; 
Ja,644) = exp! = a io “SE | (4.4) 
ca l 26 a 
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is the modified Bessel funetiom of the first kind of order & > In (4.4), ¢, Isseone 
sidered to be the strength of the direet component of s,(t) and 2¢7 is the mean- 
Squared value of the diffuse component in s,(7). For convenience deme 
“i; = oz/20- as the direct-to-diffuse power natio of s.(7). Ui 7; approachicsiaaas 
(absence of direct components). the density f, (a) represents Rayleigh fading. 
While if ys approaches infinity (presence of direct components) the density func- 
tion represents a Gaussian channel with mean «,. [Ref. 9] 

Thermal noise is assumed to be present in cach channel and is additive white 
Gaussian noise (AWGN) with zero mean and @ Warianee Of oc; = \,8. Wiciewm 
is the one-sided (AWGN) power spectral density and Bis the channel bandwidth. 
The noise is assumed to be statistically independent and identically distributed in 
each channel. For no jamming, the received Ath hop signal may be yepreseiies 


as 


Co> 


PC este On (Onan ae p(y c= 1,2 1. (4.6) 


The signal can be expanded into in-phase and quadrature components which 


gives 


r(t) = a,( cos w)t cos 0, — sin w,fsin @,) + 2,,(0) cos wt — ny (2) sin yt 
[ a, COS G. ++ Hap(t) | COS wyt = [ a, Sin 0, + nt) | Sin CO, t (4.7) 
SEPP COS (if — X47 (C) Sill Wy, 


where 71,,(¢), and n,(¢) are the independent noise quadrature components in all 
Channels. The output samples of the square-law detectors on the Ath hop are. 


dropping the time dependence, for the signal channel 


a, 2 
A ee eee 


; . : (4.8) 
= (a, COS Oy + Nex) + (ay Sin Oy + gy)”. 
and for the other channels 
2 2 
Be Geis hl = 2.3 ainee (4.9) 


When jamming occurs, it 1s considered to be present 1n all \f channels and the 
fedsiiement Channel. The jammer is assumed to spread its total available power 
eememiormly over the fraction 7 of the bandwidth 1}, O0<17<1. The average 
Femmimins Spectral density is .\,=J//' Watts Hz. The jamming spectral density 
Over the actual jammed fraction of the band is AN,/4y Watts. The combination of 
jammiing and thermal noise on the Ath hop produces the detector output samples 


for the signal channel 
> met, , 
Nik == (@;, Si] 0, a5 Meg ie) + (cd, COS Oy. ar Hex aye, : (4.10) 
and for the other channels 


Nay = (71, *K ree ar (11) ae l= ee tenis ee (4.] 1) 


uk 


where j,, and /j,,k =1,2....,£ are the independent jamming noise quadramic 
components in all channels and are assumed to be additive Gaussian noise with 
Variance defined by cs = 2/1. 

Following the derivation by Lee et al. [Ref. 5], the error probabilitvees 
pression can be derived. If equally likely M-ary symbols are assumed, the ex- 


pression for the symbol error probability, given signal of channel | is sent, is 


a 


L 
Vee) = > PPel! hops jammed) (4.12) 
/=0 
where 
L 
aE Ce — "ae (4.13) 


is the probability that /out of L channels are jammed, and 
P(e oe) (4.14) 


is the probability that an error was made given / channels are jammed. P,(c| J) is 
the probability of a correct decision given that signal of channel | is sent andy 
HODSeaTe jaumimcc. 

Whalen [Ref. 4: sec 4.7, 4.8.] shows that the square-law detector output has 
a probability density function which is noncentral chi-squared: thus, for a given 


Ain iwale i: 
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where og = NB=2N/7=2NL/F,. Whe total noise Genmin ae wes “eno. 


(Vo with probability _— n} 
Te In) wi , , (4.16) 
[ .\ a CX () aS a y/t) with probabilty i). { 


The output of the detector is then normalized by dividing by the output of the 
Maier estiiiator Which results in thewariable z, = x,/o7, i= 1,2,... Af, which 


has a conditional pdf, for the signal channel. 


6 an a (24,)'""ay, 
=n (| a ant (4.17) 
at da, ‘i 


The unconditional probability density function for z,, can be found by integrating 
ome Over all a 





Fe, (21K 1 ax) 


a5 
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or (see Appendix). 
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Where ff = 2c-R is the output signal-to-noise ratio produced by the diffuse com- 





(1 + py 


ponent. and p, = 7;A is the output signal-to-noise ratio of the specular (direct) 


cCOmponcnt of the 4 th hop signal in channel t. The ratio & is defined as 


Ve: hop energy of transmitted signal 








A=-—= : 
N total noise densitv fer 
a eee 7) 
- 7") 2:5 5) 
4 . 22 or 
The decision statistics at the output of the summers are 
L 
5= ) Fk: Pa 2 otc a Wi, (4.21) 
k=] 


Mie probability density function for =, may be found using characteristic func- 


Meee eine use O! the tables of Fourier transforms by Campbell and Foster 
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(pair 655.1) [Ref. 11: p. 79] the characteristic function of fZ, (2,,) in (4.19) is given 
by 


2 py 1 
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where p is the complex variable and 


(l + f)}-and ¢,— SS (4.23) 
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Since the z,,’S are independent and identically distributed, the characienisue 


function of the decision statistic z, can be found from 


&z,(7) =| |%z,(0). (4.24) 
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From Ref. 11: pair 650, p. 77. the probability density funcuoen for z, 1s Obtained 
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(4.26) 





@- U0. Ihe noncentral parameter of (4.26) is 
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melee: = &:/2c- 1S the direct-to-diffuse power ratio and 





is the energy per hop to noise density ratio, with £, as the symbol energy. and E£, 


@oetne Dit energy. Then. if it is assumed that /out of 2 channels are jammed. 
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Pepe tiedan direct-io-diffuse power ratio with wWW, and A; are from (4.16). 
With the use of the definition of the Bessel function (4.5). equation (4.26) can 


ipevex pressed as 
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As in the non-fading case with a fixed amplitude, the pdf for the noise only 
channels, z, i= 2.3, ... BY. is chi-squared with 2L degrees ob freedom and maggie 
obtained from (4.32) With) — 00 — a ee 


eo 


(24) = 3T(L) exp —ZiZ). 2 0) ea (4.32) 


Equation (4.14) can then be found from 
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Substituting the appropriate distributions (4.31) and (4.32) into (4.33) above 
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Equation (4.38) is a modified version of the J. P. Miller formula [Ref. 
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12: p. 42] 


for the coefficients of a composition of a formal power series and a nonunit. 


Composition is the substitution of one power series into another. 


The Miuiler 


formula applies when one series is a binomial series. Substituting (4.37) into 


(4.34) gives 
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Solving the mtegral and and using (4.41), vields 
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With some rearanging (4.42) can be 
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The portion tn the brackets 1s equal to the confluent hypergeometric Series 


fee 1; L: i) [Kef. 12: Vol I], p. 260, Ref. 14: p. 1250] and is given by! 
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] Details on the derivation of equation (4.47). 
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Li(u) 1s the generalized Laguerre polynomial, which in this case equals 
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Substituting (4.47) into (4.43) gives the desired result 
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The result from Lindsey [(41), Ref. 9] for the symbol error can be Obie 
from (4.50) by letting 7 = 0. the no-jamming case. and using the hvpergeometric 


series delinition (4.47) 
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Bit error probability from symbol error for an M-ary orthogonal svstem 
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Where the probability of symbol error is given by equation (4.12). Using (4.13), 
(4.50). and (4.52) in (4.12) gives the probability of bit error for the receiver 
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momeede non-iading case. 2¢° = 0 and «, = 1 for all A. and if the definition for the 


Laguerre polynomial (4.49) is used the result 1s equation (16) from Ref. 5 
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B. CODED FH/MFSK IN A RICIAN FADING CHANNEL WITH 
PARTIAL BAND JANIMIING 

Beicraiversii, Studied in the last Section is an effective means of overcoming 
imereetrimental effects of fading. Frequency diversity, however, is a trivial form 
Smrepetiive coding and the decision statistic of the receiver is soft decision de- 
coding of the repetitive code. More complex types of codes may also be used to 
mre ice diversity Which is much more efficient than frequency diversity [Ref. 3: 


ies |. 


4] 


Forward Error Correction (FEC) codes cam be usedaio provide impirouee 
performance with fading signals without having to use large amounts of diversity 
which cause increases in the recombining losses caused by the summing oper- 
ations. To apply coding, the circuit of Figure 8 on page 28 is modified with the 
addition of a binary to Q-ary encoder and a Q-ary to binary decoder as shown 
in Figure 10 on page 43. Information bits from a binary source at a rate RA, are 
mapped to Q-ary symbols where Q= 27 at a rate R,= R./q and applied to the 
coder input. The forward error correction coder outputs 2 Q-ary symbols for 
every A Q-ary input symbols. This gives a code rate r=A/n and the coder output 
rate is R,=a”R fk =nk,/kq. The channel modulation is FH:/MFSK with 
Af =2*. The modulator input rate is Rj, =@RjK =nk,/kKK. The M-ory eee 
are applied to a NIFSK modulator which selects one of ‘'/ baseband frequencies, 
f. i=1,2,...,.Af. The modulator output symbols are then hopped at a rate of L 
hops per svmbol, or A, = LR,=uHxLA KA. The bandwidth of the transmiuiies 
pulse is B = 2/t, where 7 = 7,/L , 7,=Ssymbol rate. The total svsten) band mares 
is made-up of ., possible frequencies spaced B Hz apart so that Ti — pea 
[Ref. 10] 

The receiver block diagram is also Shown in Figure 10 on pase 43a. 
demodulator operation is the same as described in section 4.A. After a decision 
is made as to which channel has the largest output, the M-ary words are con- 
verted into Q-ary symbols and sent to the decoder input. The decoder makes a 
hard decision as to what Q-ary svmbol was sent. A Q-ary to binary converter 
then reconstructs the binary series OB 

If the transmitted energy per information bit is defined Be (4-28), thenmume 


energy per coded M-ary word is 


AKE, ee 
Ly = 4.55) 
and the energy per hop 1s 
FE z, KKE, (4 a 
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FH/MEFSK system with coding: 


Figure 10. 


The probability of word error at the output of the decision maker is the same as 


without coding, namely 
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The noncentral parameter p, for the rate=k:n coded waveform with / out of L 
channels jammed 1s 


Rees Ep EL» : 
ar [ re +(L—J) 2 | (4.58) 


where iV; = Ay + 1\4/7 as betore. 
1. Convolutional codes 
The convolution codes considered here are binary codes with constraint 
length 7, and a rates of | 2 and 13. The Viterbi decoding aleonthimis assumes 
and the information bit error rate for the best rate, constraint length 7 code can 
be bounded by the Chernoff bound. A best rate convolutional code is one that 
cives the best performance from a Viterbi decoder. [Kel. 15] For best raicue 


COdGS ie DOUNG us 
P, <= 36D" + 211D> 2401) ees (4.59) 
and for rate 13 


Pps HD 26p ae ee (4.60) 


ad 


where 
D=2,/P. —P,). (4.61) 


In (4.61) P. is the forward transition error of a binarv symmetric channel, which 


for a binarv code and a M-arv-to-Binarv decoding is 


M/2 


2= vie ie : (4.62) 


Sete 7. iS given by (4.57). The higher order terms for D in (4.59) and (4.60) are 
neglected due to the usually small size of P. 
2. Reed-Solomon codes 

Reed-Solomon codes were found to be particarly well suited for usage on 
fading channels, more so than convolutional codes because Viterbi decoders are 
not Well-suited to the bursty errors associated with fading [Ref. 16]. The Reed- 
Solomon codes are non-binary linear block codes in which the elements of the 
code word are selected from O symbols. Usually O = 27so that & information bits 
mieeiiapped into One Osvmbol. For gq > A. where K = log,A/, q/K M-arv words 
are mapped to the Q-arv symbols and the symbol crror at the output of the de- 


coder is 
2 a Ol Saae (4.63) 


Mieoractice. g/A is usually an integer. 

For an (k.) block code. the decoder accepts an n-tuple of Q-ary coded 
svmbols and outputs a A-tuple of decoded information symbols which are further 
broken down into information bits by the Q-ary-to-binary decoder. At the de- 


coder output the probability of error for the Q-arv symbols 1s 


a Oe a3 gle Pa | 
Be mee \ ' \PAU — P.)’ ee » J Pel 2 ae (4.64) 
onl 


where dis the minimum distance between code words, which for Reed-Solomon 


codes is d=n-—k+ 1, and 


= | a | (4.65) 


is the maximum number of correctable symbol errors in a code word. [Ref. 10 ] 
The notation [x] denotes the mteger Jess than or equal to x. The probalbinaan 


bit error 1s then given by 
= (4.66) 


C. NUMERICAL RESULTS 

in both the coded and uncoded cases, the generation of the probability of bit 
error involved a computer simulation of equation (4.53). In applving optimum 
partial band jamming. the jamming fraction 7 which maximized the probability 
of bit error had to be found. This was accomplished bv performing a numerical 
search for 7 at each value of jamming power to noise poWer density [auger 
fOr Every Vale Of ik we ye andes 

1. Uncoded performance 

Uncoded performance for various levels of fading are shown in Figures 

11 to 14 for £,/.\) = 10 dB. and Af =4. Also plotted is the ideal periormiane@emaen 
no fading or jamming with the abscissa as £,/A, vice £,/\,. In COMmparine ies 
curves fur each value of # one can sce the adverse effects fading has on the System 
performance. It can also be noted how diversity improves performance signif- 
icantly for small value of 7* (Ravliegh fading). but the diversity also hinders per- 
formance when there is a strong direct component in the received signal (high 
77}, At 7+ = 15 the performance of all levels of diversity aie about the Saniewaren 
a signal-to-noise ratio of 15 dB. a similar result is seen in Figures 15 and 16. For 
the higher £,/.\, . however. the adverse effects of divcrsit aie nat realized neues 
the higher £,/.\, offsets the recombining losses. Figures 17 and I8 show a com- 


parison of performance for various values of \/ when (= 2 and £,/0\, — Genes 


46 


for both the Ravleigh fading case and a Rician fading case respectively. From 
these curves it is apparent that an increased symbol size gives better performance. 
More improvement from diversity is seen in the Rician channel. 

2. Coded performance 

The first convolutional code considered was a best rate 1/2, constraint 
length-7 code with the probability of bit error bounded by (4.59). Plots of the 
eeror rates for this code are shown in Figures 19, 20, 21, and 22 for various levels 
Gierading with Af = 4 and £,/.\, = 15 dB in worst case partial band jamming. In 
the non-fading and Rician fading with a strong direct component cases (7? = 15) 
@eeeeioures 19 and 20), the optimum level of diversity is L = 2 hops per svmbol. 
ime vicl. AS the fading becomes stronger (see Figures 21 and 22), the perform- 
ance increases when there are a greater number of hops per symbol. This is as 
expected, if these plots are compared with the results for the uncoded system in 
which the higher levels of diversity increased performance as fading increased. 
Furthermore. a comparision of Figures 20 and 22 with Figures 15 and 16 shows 
that coding does indeed decrease error rates as expected. For binary-FSK mod- 
ulation (.\/ = 2). the same gencral pattern Is evident (see Figures 23 and 2-44) where 
diversity improves performance as fading increases. Also, a comparison of the 
eee results with the MFSK results indicate that better performance is gained 
with higher values of S/. 

The best rate 1;3 constraint length-7 code results are shown in Figures 
eeeaia 26 for Ravieish fading and Rician fading when 3f = 2 and £,/N, = 15 dB. 
Again one can sce that the greater number of hops per symbol provide increased 
performance in deep fading. With the rate 1/3 code. however, slow hopping 
ff — |) provides consistently poor performance, even with no fading. A (15.8) 
meea-50lOmon code was used to generate the curves in Figures 27 and 28 for 
Af =4 and £,/:\,=15 dB in Rician and Ravleigh fading channels respectively. 
The probability of bit error for this code 1s given by (4.63) through (4.66) with 
fe orimuy=3. Therefore. the (15.8) code has an error correction capability of 
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Figure 11. Performance of receiver for no fading: y? = oo and AJ =4 under opti- 
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Figure {!4. Performance of receiver for Rayleigh fading: 3?=0 and A/=4 under 


optimum partial-band jamming conditions with number of hops-per- 
symbo! (L) as a parameter when £,/\, = 10 dB. 
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Figure 16. Performance of receiver for Rayleigh fading: ¥»?=0 and jAf=4 under 
optimum partial-band jamnung conditions with number of hops-per- 


symboi (L) as a parameter when £,/N, = 15 dB. 
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Figure 17. Performance of receiver for Rician fading: y?=15 and L=2 under 
optimum partial-band jamming conditions with number of hops-per- 


svinbol (VI) as a parameter when £,/.V, = 10 dB. 
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ficure 18. Performance of receiver for Rayleigh fading: y?=0 and L = 2 under 
optimum partial-band jamming conditions with number of hops-per- 


symbol (Mf) as a parameter when £,/\, = 10 dB. 


a 
LA 


All three codes discussed are plotted together in Figures 29 to 32 for two 
levels of diversity, L = 2 (solid line), and L = 3 (dashed line), when ‘f= 4 and 
E,/No = 15 GB. From the curves one can see that the best rate 172 convolution 
code consistently provides better detection capabilities than the other two. In 
addition, the rate 1;3 convolutional code is slightly better than the Reed-Solomon 
(15.8) code for signals with strong direct components, and much better for signals 
in deeper fading. Here again, the benifits of diversity cam be seen as the amoum 


of fading increases. 
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Fieure 19. Performance of rate 1/2 convolutional code with no fading: y? = oo and 
Af=4 under optimum partial-band jamming conditions with number 


Of heps-per- svinbol (1) as a paraineter when £,/i\, = 15 dB. 
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Figure 20. Performance of rate 1/2 convolutional code in Rician fading: 7? = 15 


and |f=4 under optimum partial-band jamming conditions with 


number of hops-per- symbol (L) as a parameter when £,/\, = 15 dB. 
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Figure 21. Performance of rate 1/2 convolutional code in Rician fading: =? = 5 and 
\f=4 under optimum partial-band jamming conditions with number 


of hops-per- symbol (L) as a parameter when £,/\, = 15 dB. 
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Figure 22. Performance of rate 1/2 convolutional code in Rayleigh fading: y? = 0 


and ‘/=4 under optimum partial-band jamming conditions with 


number of hops-per-symbol (L) as a parameter when £,/N, = 15 dB. 
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Figure 23. Performance of rate 1/2 convolutional code in Rician fading: 7? = 15 
and Af=2 under optimum partial-band jamming conditions with 


number of hops-per- svmbol (L) as a parameter when £,/.\, = 15 dB. 
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Figure 24. Performance of rate 1/2 convolutional code in Rician fading: y? = 5 and 
Af =2 under optimum partial-band jamming conditions with number 


of hops-per-svmbol (L) as a parameter when £,/N, = 15 dB. 
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Figure 25. Performance of rate 1/3 convolutional code in Rician fading: y? == 15 


and Af=4 under optimum partial-band jamming conditions with 


number of hops-per- symbol (L) as a parameter when £,/N\, = 15 dB. 
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Figure 26. Performance of rate 1/3 convolutional code in Rayleigh fading: ¥? = 0 
and ‘\/=4 under optimum partial-band jamming conditions with 


number of hops-per- symbol (L) as a parameter when £,/Nq = 15 dB. 
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Figure 27. Performance of Reed-Solomon (15,8) code in Rician fading: = y? = 15 and 
Nf =4 under optimum partial-band jamming conditions with number 


of hops-per- symbol (L) as a parameter when £,/N) = 15 dB. 
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Figure 28. Performance of Reed-Solomon (15,8) code in Rayleigh fading: y? = 0 
and Af=4 under optimum partial-band jamming conditions with 


number of hops-per- symbol (L) as a parameter when £,/A, = 15 dB. 
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Figure 29, Performance comparision of codes with no fading: 3? =oo and J/=4 
under optimum partial-band jamnung conditions with L = 2 (solid linc) 
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Figure 31. Performance comparision of codes in Rician fading: y?=5 and J/=4 
under optimum partial-band jamming conditions with L = 2 (solid line) 
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Figure 32. Performance comparision of codes in Rayleigh fading: y?=0 and 


\f=4 under optimum partial-band jamming conditions with 0= 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The theoretical study and numerical simulation of randomly spaced satellites 
indicate that there will usually be many Satellites within communication range of 
a ground terminal with relatively few periods where coverage is not available. 
The large numbers of satellites could lead to high levels of interference and many 
multipath channels if no processing of signals or directional transmission of 
signals 1s performed. Furthermore. as the range between ground terminals de- 
Ssiring to communicate increases. the number of Satellites available to communi- 
cate through decreases exponentially. This could lead to longer periods of no 
communication ability when the satellites are not equipped with a cross-link ca- 
meets. For a 250 satellite array, from five to 20 satellites can be expected 
Within range about 95°% of the time. 

The theory needed to analyze the performance of a fast-frequency-hopped 
MIFSK system with side information in a fading channel and optimum partial 
band jamming has been developed. In the configurations considered, diversity 
m@idecoding both nrovide improvement in System performance. The system 1S ex- 
pected to operate in a fading channel and will thus require diversitv and coding 
for improved performance. The optimum level of diversity for most combinations 
mieeise level and fading level is when L = 2 or L = 3 hops per symbol are used. 
Imeaddition, the best coding performance was obtained with the best rate 1:2 
convolutional code. The results for the Reed-Solomon codes mav not be an ac- 
erie vepresentation of their benefit in a bursty channel since burst analysis was 


not simulated. 


B. RECOMMENDATIONS 
Further study into other waveform configurations suitable for the LASAT 
system is needed. This study could include determining the benifits of direct se- 


@uence spread spectrum and phase shift keved (PSK) modulation. Link analvsis 


must also be performed on the proposed LASAT system. Multiuser interference, 
and a study of antenna design should be included in the analysis. The fast- 
frequency-hopping MFSK svstem theoretically developed in this thesis should 
perform well in a fading channel with jamming. To confirm the performance re- 
Sults. a hardware simulator implementing this design should be constructed. 
Network analvsis for the LASAT system should also be studied and should 


consist of the following tasks: 


e Allocating satellite capacity to voice channels in the case where voice is not 
pacxeuzea: 


© Design packet protocols to handle a mixture of packetized voice and data 
iratine: 


¢ Design highlv reliable data protocols for priority messages to be used when 
the link is forced to operate at its threshold. 


e Design connectivity algorithm for dvnamic crosslinks between LASAT’S. 
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APPENDIX 


The unconditional probability density function of z, can be derived by inte- 


grating the conditional pdf over the distribution for a, such that 


a 


Fz,,(Z1%) = | fay | A, 21k | UMA (ay)day. (A.1) 


where Fe, a,(Zin | a) and f,,(@,) are given by (4.17) and (4.4), respectively. There- 


fore, the integral to be solved is 


( 21k tH 


me, fi i) =a ! ~ f_ e | aS 
fz, (2) = 53 xP) 342 fl, Xx xP) il 7 si 502 i 


ele 
I (Z) x) Qk | QF je ve 
| ie le en icoe= 
CG 


The solution of the integral may be found by using the following indentity from 
Peer, 1/7: p. 314 














re I (anys (ond =e +? = = ) A.) 


magero J (x) is the Bessel function of the first kind. If it is noted that 


I(x) = J,x). and following substitutions are made 


po -( — a = 
20; ne (4.4) 











mie resull 1s 


2 
Y je 


Z\21k) = =n P| ee Sy ian 
wi 21 + 267/267) 2(1 + 207/2e%) JL (1 + 267/263) 


This equation can be put into a more convienent format if the substitutions of 


Px = ORR = af/2e7 and B = 26°R = 207/2c} are made 


& 


er Sime Gi) ae: 
Feu) = FFB) 0) (1 + fi) | ar (4.6) 


where p, is the output signal-to-noise ratio of the specular (direct) component, 


and f is the output signal-io- noise ratio of the diffuse component. 
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